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use). A second system provides several aviation applications 
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craft separation, flight arrival management, flight re-routing, 
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AIR TRAFFIC MANAGEMENT EVALUATION 
TOOL 

ORIGIN OF THE INVENTION 

The invention described herein was made, in part, by one or 
more employees of the United States Government and may be 
manufactured and used by or for the Government for govern- 
mental purposes without the payment of any royalties thereon 
or therefor. 

TECHNICAL FIELD 

The present invention is a method and system for evaluat- 
ing and implementing selected air traffic management con- 
cepts and tools. 

BACKGROUND OF THE INVENTION 

In the United States, as many as 7,000 commercial and 
private aircraft may be in the air simultaneously at a given 
time and date, and the total number of commercial flights in a 
given 24 -hour period generally exceeds 50,000. For example, 
in March 2001, more than 57,000 fli gilts were reported for 
one 24-hour period. Further, the growth in commercial air- 
craft traffic has been growing at a rate of between 2 and 7 
percent per annum. Faced with a doubling of commercial air 
traffic in a time interval of between 10 and 35 years, workers 
in aviation are concerned with implementing air traffic man- 
agement approaches that can safely and reliably handle air 
traffic growth over the next several decades. 

What is needed is an approach that receives proposed flight 
plans and associated flight route information and flight 
parameters for a plurality of aircraft operating in a given 
region (e.g., the continental United States) and provides 
actual flight routes and schedules, based upon expected air 
traffic, and that avoids or minimizes air traffic incidents, by 
changing one or more flight plan parameters where appropri- 
ate, for one or more of these aircraft. Preferably, the system 
should provide flight route information and parameters for 
normal flights, for direct-to flights, for emergency responses 
and for free flight responses to events. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
method and system for evaluating and implementing air traf- 
fic management (ATM) tools and approaches for managing 
and for avoiding an air traffic incident enroute, before the 
incident occurs. The invention includes a first system that 
receives parameters for flight plan configurations (e.g., initial 
fuel carried, flight route, flight route segments followed, flight 
altitude for a given flight route segment, aircraft velocity for 
each flight route segment, flight route ascent rate, flight route 
descent route, flight departure site, flight departure time, 
flight arrival time, flight destination site and/or alternate fight 
destination site), flight plan schedule, expected weather along 
each flight route segment, aircraft specifics, airspace (alti- 
tude) bounds for each flight route segment, and navigational 
aids available. The invention provides flight plan routing, 
direct routing and/or wind-optimal routing, using great circle 
navigation using spherical Earth geometry. The invention 
provides for aircraft dynamics effects, such as wind effects at 
each altitude, altitude changes, airspeed changes and aircraft 
turns to provide predictions of aircraft trajectory (and, option- 
ally, aircraft fuel use). 
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A second system provides several aviation applications 
using the first system. Several classes of potential incidents 
are analyzed and averted, by appropriate change enroute of 
one or more parameters in the flight plan configuration, as 
5 provided by a conflict detection and resolution module and/or 
traffic flow management modules. These applications include 
conflict detection and resolution, miles-in trail or minutes-in- 
trail aircraft separation, flight arrival management, flight re- 
routing, and weather prediction and analysis. 

In one approach, the present flight plan configurations for 
each of two or more aircraft are analyzed, and the system 
determines if an aircraft flight conflict (distance of closest 
approach of two aircraft less than a threshold number, such as 
15 3-8 nautical miles) is likely to occur during or at the end of the 
flight of the aircraft. If occurrence of a conflict is likely, the 
system remodels the flight plan configuration(s) for one or 
more of these aircraft, analyzes the remodeled 
configuration(s), and determines if a conflict is likely with the 
20 remodeled flight plan configuration(s). If the answer to the 
query is “no,” the system accepts and optionally implements 
the remodeled flight plan configuration(s) for the aircraft 
flights being examined. If the answer to the query is “yes,” the 
system further changes one or more parameters in the remod- 
25 eled flight plan configuration(s) and again inquires if a con- 
flict is likely to occur with the changed and remodeled flight 
plan configuration(s). This procedure is iterated upon until a 
remodeled flight plan configuration is found that avoids a 
conflict along the flight route. Changes to be made to avoid a 
30 conflict may be split between the two aircraft, or allocated to 
a single aircraft, according to a selected sharing fraction § 
(O^^l). 

In another approach, the system analyzes consecutive air- 
craft spacing along a selected flight route segment. If the 
35 spacing for two consecutive aircraft is smaller than a thresh- 
old number, the relative velocity of one or both of the two 
aircraft is adjusted to maintain at least the threshold spacing. 

In another approach, the system analyzes flight arrival 
40 information for a selected destination (airport) and deter- 
mines if the destination will be too congested when a selected 
aircraft arrives there at its scheduled arrival time. If the 
answer to the query is “yes,” departure of the selected aircraft 
is delayed by an appropriate time interval so that an arrival 
45 slot for the aircraft is likely to be available at the now-modi- 
fied estimated time of arrival. 

In another approach, the system analyzes weather informa- 
tion along a selected flight route to a selected destination 
(airport) and determines if the anticipated weather is too 
50 severe. If the weather along the selected flight route is too 
severe, (1 ) the remainder of the flight route is altered to arrive 
at the same destination or (2) the remainder of the flight route 
is altered to arrive at an alternative destination. Flight route 
alteration can be implemented enroute or before departure. 
55 The system relies upon several integrated and interacting 
modules. In a first module, a flight route is specified, as a 
sequence of waypoint locations and altitudes or as a route 
specified in the National Playbook Routes or in the Coded 
Departure Routes. In a second module, flight route and air 
60 speed restrictions are imposed, as determined from a miles- 
in-trail or minutes -in-trail restriction (“MIT” restriction), a 
ground delay restriction and/or a ground stop restriction. A 
third module provides individual aircraft rerouting around a 
congested area and a fourth module to avoid a conflict with 
65 another aircraft, in which the predicted nearest distance of 
approach of the two aircraft is less than a selected threshold 
distance. 
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The core system can be operated in at least five modes: (1) 
a playback mode, in which stored data from earlier flights or 
runs is played back for evaluation and further analysis; (2) a 
trial planning mode, in which selected parameters are altered 
and one or more situations are re-run to evaluate the impact of 5 
these alterations; (3) a simulation mode, in which filed flight 
plans and modifiable initial conditions are used to predict 
aircraft locations and to forecast or predict traffic patterns as 
a function of time; (4) a live mode, using filed flight plan and 
tracking information collected by air traffic controllers to to 
provide aircraft locations in real time; and (5) a batch or 
collective mode, to provide a consolidated view or probabi- 
listic view of the collective effects of variations in several 
initial conditions, parameters and scenarios. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates architecture of a server according to the 
invention. 

FIG. 2 illustrates components of a core architecture accord- 20 
ing to the invention. 

FIG. 3 illustrates a three dimensional screen display of 
NAS flights enroute, indicating ascent of each flight. 

FIG. 4 illustrates effect of local wind on aircraft heading. 

FIG. 5 illustrates a GUI screen, according to the invention, 25 
displaying NAS flights enroute within the continental con- 
tiguous U.S. at a particular time. 

FIG. 6 illustrates geometrical and physical parameters of 
concern in an aircraft flight. 

FIG. 7 illustrates two aircraft traveling along the same 30 
route segment. 

FIG. 8 illustrates two aircraft traveling in the same region. 

FIG. 9 illustrates a conflict situation for two aircraft. 

FIG. 10 illustrates direct-to routing. 

FIG. 11 is an example of a display of National Playbook 35 
Routes between major airports on the West Coast and on the 
East Coast. 

FIG. 12 illustrates rerouting of east-bound and west-bound 
flights around a convective weather cell. 

FIG. 13 graphically illustrates cumulative aircraft delay 40 
contours resulting from joint time delays in departure rates 
from two adjacent airports. 

FIGS. 14 <3, 14 b, 14 c and 14 d schematically illustrate an 
embodiment of a procedure for practicing the invention. 

45 

DESCRIPTION OF APPLICATIONS OF THE 
INVENTION 

FIG. 1 illustrates the architecture of the system, emphasiz- 
ing sources of the information used by the system. A geo- 50 
graphically distributed or central server group 11 includes a 
route parser and trajectory modeler module 13, an air traffic 
analyzer module 15 and a graphical user interface (GUI) 17. 
The server group: receives weather information from the 
National Oceanics and Atmospheric Administration 55 
(N.O.A.A.) and/or from the U.S. Weather Bureau 21; receives 
aircraft flight path and location information from the F. A. A. ’ s 
enhanced traffic management system (ETMS) 23; receives 
aircraft performance data, including aircraft climb, cruise and 
descent information, from an aircraft performance database 60 
25; and receives flight adaptation information on airports, 
airways, and traffic control centers and sectors from a flight 
adaptation module 27. The server group 11 analyzes the 
received information and provides at least six types of out- 
puts: (i) flight deck-based conflict detection and resolution 65 
(CD&R); (ii) airport arrival and departure rules (iii) direct-to 
routing analysis for use in planning direct-to flights; (iv) air 
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traffic integration information; (v) evaluation of an initial 
playbook route and subsequent changes that have been or will 
be implemented; and (vi) system-wide optimization of flight 
routing, flight departures and flight arrivals. The system 
focuses upon flights for which a flight plan has been filed 
(referred to as “NAS flights” herein). The system relies upon 
a combination of: (1) several relevant and periodically 
updated databases that provide information on aircraft con- 
figurations and performance data, locations and configura- 
tions of available airports and runways, special use or 
restricted airspaces, and present and estimated fiiture weather 
data; (2) software applications that provide computations, 
forecasting and/or visual presentations; (3) a GUI that pro- 
vides static and/or animated views of present and/or predicted 
air traffic, in a selected airspace region, Air Route Traffic 
Control Center (ARTCC), ARTCC sector and/or nationwide; 
and (4) an output signal stream providing recommended con- 
trol advisories for traffic flow specialists. 

In one embodiment, the GUI 17 provides: (1) an option of 
two dimensional or three dimensional displays of a particular 
aircraft configuration in a region; (2) separate or integrated 
displays of air traffic, wind components, weather and/or adap- 
tation elements; (3) animated displays of three dimensional, 
weather and/or air traffic forecasts; (4) displays of filtered air 
traffic as presented, using traffic stream visualization to sup- 
press display of selected classes of air traffic; and (5) fly-by 
animated displays, using a scroll bar to view past, present and 
future positions and conditions of air traffic and weather 
patterns. 

FIG. 2 illustrates the architecture of the core components of 
a route parser and trajectory prediction module 13 for the 
system. This module provides wind data 31 and information 
from a route navigation module 33 to determine aircraft head- 
ing commands, which are received by a heading dynamics 
module 41. The heading dynamics module optionally 
includes information on maximum banking angle at one or 
more altitudes and maximum turn rate at one or more alti- 
tudes. The route navigation module 33 receives information 
from a direct routing module 35 or, alternatively, from a flight 
plan routing module 37 and provides destination coordinates. 
An airspace module 39 provides information to a flight option 
logic module 40 that determines whether the flight is simu- 
lated according to direct-to routing or according to flight plan 
routing. Where a flight plan is filed and followed, the flight 
plan routing module 37 may provide coordinates of one or 
waypoints for the flight route. 

An aircraft performance database 44 provides relevant per- 
formance information on more than 500 aircraft, optionally 
including data for each aircraft on maximum airspeed in 
absence of wind, fuel consumption at different altitudes, dif- 
ferent air speeds and different payload weights, maximum 
climb rate at one or more altitudes, aircraft weight range 
(empty to fully loaded), practical maximum flight altitude, 
and angle of attack at initiation of stall (optional). This infor- 
mation is provided for and used by an aircraft performance 
module 45 that models a selected aircraft’s performance and, 
in turn, provides airspeed command and performance limits 
information for an airspeed dynamics module 47. The aircraft 
performance module 45 also provides altitude command and 
performance limits information for an altitude kinematics 
module 49. The airspeed dynamics module 47 provides rel- 
evant, processed airspeed and altitude information to the lati- 
tude and longitude kinematics command module 43 and to 
the heading dynamics module 41. The latitude and longitude 
(LLK) module 43 also receives relevant, processed informa- 
tion from the altitude kinematics module 49 and information 
on flight path angle. The wind data module 31, the airspace 
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module 39, the aircraft performance module 45, the LLK 
module 43 provide output information that is received by the 
graphical user interface 17. 

A. Graphical User Interface (GUI) 

The GUI 17 optionally provides a three-dimensional view 
of one or more selected ARTCC sectors, an ARTCC itself, a 
geographic region, or the continental contiguous U.S. or 
Alaska or Hawaii, as illustrated in FIG. 3, in which the view 
is from the side, not the top, and an aircraft climb path or 
descent path is represented by an almo st-vertical line in this 
view. 

The GUI 17 can display winds-aloft patterns at selected 
altitudes (e.g, FL180, FL 230, FL 270, FL 310, FL 350, FL 
410 and FL 450), corresponding to well-used cruise altitudes 
for commercial flights, for one or more selected ARTCC 
sectors, an ARTCC itself, a geographic region, or the conti- 
nental contiguous U.S. or Alaska or Hawaii. The GUI can also 
display weather patterns, horizontally and vertically, which 
have developed or are likely to develop along a selected flight 
route or in a sector or an ARTCC, optionally using color 
coding or texture coding to display different adverse or 
unusual weather conditions. 

The three dimensional, weather and NAS air traffic forecast 
visual presentations can be animated for update and display at 
time intervals of 1 -60 minutes. The air traffic stream can be 
filtered so that only a relevant portion of the NAS air traffic is 
displayed, or is displayed in a different color or other indi- 
cium, based upon parameters such as airline (commercial 
flights only), aircraft manufacturer, aircraft capacity, flights 
within a selected heading angular sector, flights within a 
selected altitude band, flights having a selected source, flights 
having a selected destination, or flights having an estimated 
time of arrival (ETA) within a selected time interval at a 
selected destination or group of destinations. This filtering 
capability is useful for estimating or visualizing the airport 
arrival demand at a selected destination and for visualizing 
enroute flight segment and airport demand, within a specified 
time interval. 

B. Provision and Evaluation of Weather and Winds Data 

Assessment of weather date (including winds) at various 

altitudes is integrated into the system, using weather and/or 
wind information sources such as Collaborative Convective 
Forecast Product (CCFP), NOWRAD, National Convective 
Weather Forecast (NCWF) and Corridor Integrated Weather 
System (CIWS). CCFP and NCWF are national scale weather 
forecast products that are provided by the Aviation Weather 
Center. CCFP provides two-hour, four-hour and six-hour 
forecasts that are updated every two hours, and NCWF pro- 
vides an hourly forecast. CIWS is a high resolution weather 
forecasting product that focuses on the northeast region of the 
United States and provides storm location information, echo 
tops and an animated two-hour forecast for growth and decay 
of storms. NOWRAD, developed by Weather Services Inter- 
national, provides high quality national and regional radar 
imagery. The system also allows a user to identify flights that 
are projected to fly through one or more specified CCFP- 
defined weather cells and to automatically provide a re-rout- 
ing for selected flights that are adversely impacted by weather 
in such cells. A Rapid Update Cycle (RUC) winds module, a 
product of the N.O.A.A., is used in the trajectory prediction 
module of the system, and a wind-optimal re-routing algo- 
rithm is utilized to estimate the most ftiel-efficient route(s) 
between a source and a destination. 

Optionally, the system provides optimal routing in the 
presence of wind and/or flight constraints. In a relatively 
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uncomplicated embodiment, for a single leg or segment in a 
flight route, if the local wind at the anticipated cruise altitude 
has a velocity vector v w =(v w cos 0 w ,v w sin 0 W ) and the aircraft 
has a true air speed of v^ and is to travel at an angle Q a , comp , 
5 relative to true north or magnetic north, after accounting for 
the effects of wind, the thrust of the aircraft should be oriented 
at a modified angle Q a>comp , given by 

tan 0hd a,comp=(sin Q a , comp -p sin 6 w )/(cos d apomp - 
10 p cos 0 W ), (1) 

P =yjy a , (2) 

as illustrated in FIG. 4. The aircraft true air speed is estimated 
15 by 

Va={v«, C0W /+V w 2 +2v a ^ 0 ^V w COS(0^ o „^-0J} 1/2 (3) 

C. Interpolation of Wind and Weather Data 
20 Each weather variable (including wind variables), collec- 
tively denoted W(x, y, z, t), is measured at a relatively small 
number of spaced apart locations and at times that are sepa- 
rated by one to six hours or more. An aircraft flight crew will 
25 need to estimate a value of the variable W at a location that is 
spaced apart from the measurement location and at a time that 
does not coincide with any measurement times for that vari- 
able. The system optionally provides an estimation procedure 
that interpolates between the measured values at the measure- 
30 ment locations to provide a continuously varying function 
value that coincides with each of the measured values at the 
measurement locations. Let {r„}„ be a sequence of spaced 
apart location vectors corresponding to the measurement 
locations, r w =(x w , y w , z w ) for the variable W(r,t) at the most 
35 recent time(s) the variable W was measured. Each set of four 
nearest neighbor location vectors {r w } w defines a tetrahedron, 
having the location vectors as vertices, and the collective set 
of tetrahedrons fills all space, with overlap at boundary planes 
40 for any two contiguous tetrahedrons. 

Ignore the time variable t and consider a location vector 
r=(x,y,z) lying in the interior or on a boundary of a selected 
tetrahedron Te(l, 2, 3, 4) defined by four spaced apart, non- 
coplanar measurement location vectors, r w =(x M , y w , z n ) (n=l, 
45 2,3,4), at which the measurement values W (r„)=W (x„ , y„, z„) 
are known. The estimation function 

W(r, est ) = ( 4A ) 

WftiHk - r 2 \\r - r 3 \\r - r 4 |}/{ki -r 2 |ki -r 3 |ki -r 4 \} + 
^(r 2 )q|r-r 1 ||r-r 3 |k-r4|}/lk 2 -r 1 |k 2 -r3lk 2 -r 4 |}4- 
W(r 3 ) -\\r-nWr- r 2 ||r - r 4 \}/{\r 3 - n\\r 3 - r 2 \\r 3 - r 4 1} + 

55 W(r 4 ) -Ik - Hlk - r 2 \\r - r 3 \}/{\r 4 - n ||r 4 - n\\?A ~ *>ll 


is continuous within the tetrahedron Te(l , 2, 3, 4) and satisfies 
60 W(r=r„;est)=W(r w ). Because the measurement locations are 
spaced apart (in at least one of the three coordinates x, y and 
z), the denominators in Eq. (4) are never 0, and the magnitude 
of the function W(r;est) is bounded. The enveloping figure 
Te(l , 2, 3, 4) can be extended to a general polyhedron, includ- 
65 ing a line segment, a triangle, a tetrahedron and any polyhe- 
dron having two or more boundary surfaces (endpoints or 
vertices). More generally, if measured values W(r„) are pro- 



US 7,702,427 B1 


7 

vided at N distinct points, r=r„ (n=l, . . . , N; N^4), a suitable 
estimation function is 


W(r; est ) = \ W(r m ) Y\ “ r R )l(r m - r„)}. 

‘ 1 n=l,n+m 


Where the location vector r lies within or on a triangle Tr(l , 
2, 3) defined by three spaced apart, non-collinear measure- 
ment location vectors r’ w . (n ,= l , 2, 3) that serve as vertices for 
the triangle, the estimation function may be expressed as 

w'{n est ) = -rfiki -rfl} + (5) 

w(/ 3 ) ■ Ik - dik - 'ill/M - dlks - f 2 |}, 

where the interpretations are similar to those for the estima- 
tion function W(r;est) in Eq. (4). 

Where the location vector r lies on a line segment Ls(l, 2) 
defined by two spaced apart measurement location vectors 
r" w . (n"=l , 2) that serve as endpoints for the line segment, the 
estimation function may be expressed as 

W" (r; est) = (6) 

W(r'{) dk - r$\ }/K - >2 II + Wtf) ■ {\r - rfll/K - Z[\}, 


where the interpretations are similar to those for the estima- 
tion functions W(r;est) and/or W'(r;est) in Eqs. (4) and (5). 

More generally, one can define an estimation function 
W* (r;est) as a sum of two or more continuous characteristic 
functions W*(r;k) (k=l, . . . , K'; K'i^2), where the character- 
istic function W*(r;k) satisfies 


W*(r=r p -k)=W(r k ) (p = k) (7A 

= 0 (P * k). (7B 


The function W(r;est) or the function W*(r;est) allows inter- 
polation of a weather-wind value for any location within a 
polyhedron of dimension 1 or higher, defined by measure- 
ment location vectors as vertices of the polyhedron. 

The values W(r„) in Eq. (4) can be replaced by time- 
dependent weighting functions W(r„;t-t„) that are monotoni- 
cally decreasing with the time difference, t-t„, (i^O) between 
the present time t and the (most recent) time t„ at which the 
measurement W(r„) was taken. An example of such weighting 
functions is 

W(r n -t~t n )= (8) 

Pn ■ W{r n ) ■ exp {-a n {t - r„)} + (1 - fi n ) ■ W(avg){l - exp{-ar„(r - r„)}} 

where a n is a small positive first selected weighting index, 
is a second selected weighting index satisfying 1 , and 
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W(avg) is a suitable representative value of the variable W for 
a location associated with the vector location r. 

D. Wind Optimal Routing and Other Route Choices 

5 A system user can choose among any of three or more 
routing procedures: (1) a user-preferred route between two 
waypoints, including but not limited to a route from origin 
airport to destination airport; (2) an NPR Direct route, which 
uses a National Playbook Route; and (3) a wind optimal route, 
to as disclosed in U.S. Pat. No. 6,600,991, J incorporated by 
reference herein. In one embodiment, a “wind optimal route” 
is determined by (i) providing a nominal route between first 
and second waypoints in the presence of a first wind environ- 
ment; (ii) providing values for a second wind environment 
15 that differs from the first wind environment; and (iii) using a 
computer to determine a neighboring optimal control solution 
for an aircraft moving at a selected speed between the first and 
second waypoints in the presence of the second wind envi- 
ronment. In one approach, the neighboring optimal solution 
20 provides a differential solution that determines one or more 
route increments that suffice to move the aircraft from the first 
to the second waypoint when the first wind environment is 
modified to become the second wind environment. The dif- 
ferential solution may be expressed in terms of latitude and 
25 longitude coordinates, in terms of modifications to a great 
circle route, or in other terms. 

E. Use of Filed Flight Plans 

The system receives and stores a flight plan for each NAS 
30 flight, which includes all flights governed by instrument flight 
rules (IFR), for which a flight plan must be or is filed. Flights 
for which a flight plan is not filed are not covered by the 
system. The GUI 17 , working in combination with other 
modules, provides a two-dimensional top view of NAS air 
35 traffic, with each aircraft being represented by a visually 
perceptible symbol, such as a cross or a generic plan view of 
an airplane. Optionally, different types of aircraft can be 
represented by visually distinguishable symbols (e.g., in dif- 
ferent colors, different sizes or different symbols; commer- 
40 cial flights versus other NAS flights). The NAS air traffic can 
be illustrated for one or more selected sectors of an ARTCC 
(22 at present), an ARTCC itself, a geographic region, or the 
continental contiguous U.S. or Alaska or Hawaii. Each 
ARTCC may have each staffed by a team of air traffic con- 
45 toilers (ATCs). FIG. 5 illustrates a GUI screen showing 
approximately 4530 aircraft enroute within the contiguous 
states at a particular date and time (18 Mar. 2000 at 20:26 
UCT). The system can provide views similar to FIG. 5 at time 
intervals of 1-60 minutes, or longer if desired, using aircraft 
50 location predictions determined from the flight plan. 

When a flight plan is altered by the appropriate ATC, the 
flight plan alteration will normally be electronically posted to 
the ETMS and will be picked up by the system. The extant 
flight plan is then altered accordingly in the system flight plan 
55 database. 

F. Aircraft Performance Database 

Aircraft performance parameters for more than 500 repre- 
sentative aircraft models are provided in an aircraft perfor- 
60 mance database, currently provided by the Base of Aircraft 
Data (BADA), developed and maintained by the Euro Central 
Experimental Center in France, which is part of the system. 
Table 1 illustrates the parameters available for a representa- 
tive aircraft, a Boeing B757. The Table first provides cali- 
65 brated air speed schedule for a standard CAS -Mach climb 
(290 knots calibrated air speed to Mach 0.78), for a standard 
cruise rate (320 knots or Mach 0.80) and for a standard 
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descent rate (300 knots CAS or Mach 0.78). As altitude 
increases, the true air speed (TAS) increases faster than indi- 
cated air speed (IAS). 

Table 1 also sets forth cruise data for different flight levels 
FL=30-420 (MSL altitudes of 3,000-42,000 feet), a corre- 
sponding optimum TAS for that FL, and fuel consumption 
(kgm/min) for each of three aircraft mass loading configura- 
tions, m=69,600 kgm (low mass), m=95,000 kgm (nominal or 
medium mass) and m=l 10,000 kgm (high mass). TAS 
increases monotonically with altitude or flight level to a cer- 
tain Mach number, then decreases and subsequently levels off 
with further increases in altitude. Fuel consumption varies 
markedly with altitude, especially for a high mass configura- 
tion. 

Table 1 also sets forth optimal climb or ascent rate at flight 
levels FL=0-420 for low, medium and high mass loading 
configurations. Table 1 sets forth optimal descent rates at 
flight levels FL=0-420, for a medium mass loading configu- 
ration. Table 1 is an example of the aircraft-performance data 
for more than 500 aircraft that are included in the system. 

The ascent rates and descent rates set forth in Table 1 are 
recommended rates for all altitudes. For altitudes above the 
transition altitude (normally between 15,000 and 20,000 feet 
MSL), the ascending or descending aircraft may follow a 
programmed altitude rate change. 

An aircraft ascending to a cruise altitude will often follow 
one of a set of specified programs of air speed and climb rate. 
The programs may include a prescription for maximum climb 
rate (referred to as V x ) and/or a prescription for maximum 
angle of climb (referred to as V ), as well as other special 
purpose ascent rate prescriptions. 

An aircraft making a constant rate turn will have a turn rate 
limited by the allowable stress, the aircraft air speed, the 
density altitude and other relevant variables. Turn rates are 
typically in a range of 1 -4 degrees/ sec For example, a turn rate 
of co=3 degrees/sec (0.05236 radians/sec) requires 120 sec to 
execute a 360° turn. 

G. Airports, ARTCC Sectors and Air Traffic Monitoring 

The system applies NAS air traffic demand forecasting and 
management to provide flight planning and/or replanning, for 
example, through change of destination, change of cruise 
altitude, change of cruise speed or change of flight 
waypoint(s), to comply with an applicable MIT flight restric- 
tion or a flight separation requirement that is implemented. 
This may include restrictions based upon airspace class and/ 
or special use airspaces. The system provides on-demand 
reports of number of NAS flights that are known to be within, 
or are predicted to be within, a specified ARTCC, an ARTCC 
sector, a flow constrained area (FCA) and/or a special use 
airspace (SUA), at a selected time or within a selected time 
interval, using historic, stochastic, forecast and/or determin- 
istic models of the NAS flights. Presently, 22 ARTCCs and 
about 830 ARTCC sectors are defined, and a given ARTCC 
may have a super-high (altitude) sector overlying one or more 
high sectors and a high sector overlying one or more low 
sectors. 

The system can be used to design efficient aircraft ground 
delays and/or ground stops at a selected airport. The available 
visual displays include screen displays, histograms, bar 
charts, tables and map displays. 

Where an ARTCC sector or a special use airspace (SUA) or 
a flow constrained airspace (FCA) experiences increased or 
unusual demand, this sector or SUA and adjacent regions may 
be rearranged or reformatted, for example, (i) by decompos- 
ing the affected sector or SUA into two or more sub-regions, 
each with its own air traffic controller (ATC) or set of flight 
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restrictions and/or (ii) by rearranging the boundaries of the 
region and adjacent regions to balance the load on the ATC 
assigned to each of the regions. The system allows manual, 
visual modification of ARTCC sector boundaries and special 
5 use airspace boundaries and integrated display of air traffic 
within these modified boundaries. Modified and unmodified 
boundaries and air traffic can be displayed in two and three 
dimensions, with optional playback, simulation and live pre- 
sentations. Sector, SCA and FCA demand reporting can be 
to visualized using this option. Using any of the available sys- 
tem display modes (live, playback or simulation), display of 
NAS air traffic through the sector or SUA or FCA can be 
manually modified, using an intuitive click-and-drag capabil- 
ity built into the GUI component to implement a what-if 
1 5 scenario that displays the results of reconfiguration of a sector 
or an SUA. Two dimensional and three dimensional visual- 
izations and air traffic reporting are available for the 
(changed) sector and/or SUA and/or FCA boundaries and for 
the resulting (re)allocation of air traffic. The predicted 
20 demand on thus-modified NAS resources can thus be mod- 
eled and analyzed, using selected air traffic flow metrics. 

H. Route Parser and Trajectory Predictor 

FIG. 6 illustrates some geometric and physical parameters 
for an aircraft in flight. The aircraft has a present location 
vector 

r=(r cos X-cos t,rcos X-sin T,rsin X) (9) 

and moves with a present velocity vector (ignoring wind 
30 effects) 

v=(vcos a-cos p,vcos a- sin fJ,vsin a), (10) 

where r and v are the aircraft radius vector and velocity vector, 
measured relative to the Earth’s center. Here, t and X are 
35 longitudinal and latitudinal angles, respectively, measured 
from a reference position, such as the prime meridian and/or 
the equatorial line, and a and p are velocity vector angles. 

An LLK module in the invention utilizes spherical Earth 
equations of motion for an aircraft, 


3X/3f={v cos X-cos t+w^/R, 

(11) 

3t/3?={v cos X-sin t+w e }!{R cos X), 

(12) 

sin -1 {(8/z/3r)/v}. 

(13) 

r(X,T;r)=r(Earth;mean)+/z (X,T;f) , 

(14) 


where w^ and w E are the north-directed and east-directed 
components of local wind velocity, x is longitudinal or azi- 
50 muthal angle for the aircraft location, X is latitude or polar 
angle for the aircraft location, and h=h(L,x;t) is AGL height 
(measured relative to local ground level, rather than relative to 
sea level) of the aircraft above the local terrain. 

Using the system, creation of portions of air traffic sce- 
55 narios can be automated, partly relieving an air traffic mod- 
eler of what would otherwise be a manually intensive proce- 
dure. Filtering and historical flight plan databases associated 
with the system can be used to extract historical air traffic 
patterns (optionally, over two or more flight days) from 
60 archived data, for flight plans that were followed and for 
deviated flight plans. An intuitive flight creation GUI allows 
flights to be added to (or deleted from) the historical air traffic 
patterns. The scenario creation module can be used to develop 
futuristic air traffic scenarios that will conserve scarce NAS 
65 resources. 

Optionally, certain of the computations and the displays 
can be abbreviated or simplified in order to allow NAS flight 



US 7,702,427 B1 


11 

modeling on a laptop computer, using a parametric trajectory 
prediction engine, as opposed to modeling on a more elabo- 
rate (and less portable) computer system. A simplified flight 
trajectory prediction model may use linear trajectory predic- 
tion or may use a more elaborate quadratic trajectory predic- 
tion, in which a great circle route is approximated, as dis- 
cussed in Section K. The system architecture uses a 
combination of Java and C coding and can work in the Macin- 
tosh, Windows, UNIX and LINUX platforms. 

I. Traffic Analyzer 

The system enables demand forecasting of air and ground 
traffic to predict or estimate (1 ) number of flights in a selected 
sector, (2) number of flights along a selected segment of a 
flight route or airway, (3) airport arrival and departure rates, 
(4) demand for selected special use airspaces and (5) demand 
for flow constrained areas. 

A fleet impact assessment module allows a user to deter- 
mine if a selected flight in an airline’s schedule will be 
impacted by a specified NAS constraint. The constraint may 
be a weather cell, an active special use air space, a congested 
resource (e.g., a sector, an airway, an airport or a particular 
runway. A special display screen optionally displays the 
impacted flight, relevant details of the associated flight plan 
and the NAS constraint. Optionally, a potential impact of the 
constraint on an alternative flight plan can also be demon- 
strated. 

The system provides demand forecasting concerning the 
number of flights, airports, sectors, special use airspaces and 
flow constrained areas. Demand is predicted based on a com- 
bination of stochastic modeling, forecasting, deterministic 
modeling and/or actual historical counts and can be coupled 
with models of traffic flow management restrictions or con- 
straints (re-routing, ground delay, ground stop, and miles-in- 
trail and minutes -in- trail (“MIT”) restrictions. Displays of 
forecast variables are available as bar charts, tables and map 
displays. 

If a landing slot is likely to be available for the selected time 
interval at the selected destination, the system advises that the 
flight can proceed as planned. If a landing slot is not likely to 
be available in the selected time interval at the selected des- 
tination, or if the weather along at least a portion of the 
planned flight route is likely to be too severe, the system 
advises the aircraft of the slot non-availability and/or inclem- 
ent weather and optionally: (1) provides an alternate destina- 
tion for the flight where a landing slot will be available during 
a corresponding time interval of arrival (“TIOA”); (2) advises 
delay of departure of the flight until a time corresponding to a 
time-delayed TIOA, when a landing slot will be available; (3) 
selects an alternative destination (for the enroute aircraft), 
consistent with the remaining fuel reserve for the aircraft and 
existing weather along the alternate route, for which a landing 
slot will be available at a corresponding TIOA; and/or (4) 
advises postponement or cancellation of the flight. The sys- 
tem optionally estimates the remaining fuel for the aircraft, 
before directing the aircraft to an alternative destination. 

J. Miles -in-Trail and Minutes-in-Trail Restrictions 

FIG. 7 illustrates a spatial relationship between first and 
second aircraft (n=l and n=2) traveling consecutively along 
the same route segment RS . The two aircraft need not have the 
same departure site or the same destination site. All that is 
required is that the two aircraft travel the same route segment 
for a portion of the total route of each aircraft, within a given 
time interval having a time interval length, such as At(seg- 
ment)=2-7 min. According to an MIT restriction, the two 
consecutive aircraft are required to maintain either (1) a mini- 
mum distance of separation d(thr)=3-50 miles along the route 
segment (miles-in-trail), depending upon the present loca- 
tions of the two aircraft, or (2) a minimum temporal separa- 
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tion At(thr), typically 0.6-3.33 minutes (minutes-in-trail). For 
a given initial time t=tl, an initial location vector r M and an 
initial velocity vector v l z - is determined for each of the air- 
crafts, i=l, 2. A separation distance along the common route 
5 segment 

d{t)=\r Ui +v Ui {t-t\)-r U2 -v U2 {t-tl)\ (15) 

is then determined, using a linear approximation, for all times 
{tl^t^t(sep)} for which both aircraft will remain on the 
io common route segment, where the vectors v 1;1 and v 12 are 
parallel but do not necessarily have the same magnitude. The 
calculation of minimum separation distance, given by 

^(min) 2 =Ar l52 2 Av 1?2 2 -(Ar u -Av l52 ) 2 }/(Av 1?2 ) 2 , (16) 

15 and the calculation of time of minimum separation distance 

/(min)-n=-(Ar 1?2 -Av lj2 )/(Av 1>2 ) 2 , (17) 

are analogous to those for the FIG. 2 configuration but is more 
straightforward because v 1?1 and v 12 are parallel in this situ- 
20 ation. If d(min)sid(thr) and 0^it-tl ^t(sep)-tl , the system 
notifies one or both aircraft and requests that at least one of the 
two aircraft change at least one of the parameters of the 
velocity vector(s) v ljZ - (i=l, 2). If, for example, aircraft no. 1 
precedes aircraft no. 2 and v x 1 *v 1 2 <v i 2 ' v i 2 ? (1) the second 
25 aircraft can reduce its speed lv 12 l, (2) the first aircraft can 
increase its speed W 1 x I, (3) one of the two aircraft can change 
its flight altitude (usually, by a multiple of 2000 feet), or (4) 
one of the two aircraft can change its flight route, and (5) one 
of the two aircraft can change its flight departure time (if at 
30 least one of the two aircraft has not yet departed) so that the 
separation distance d(t) does not decrease to or below d(thr) 
during the time interval {tl^iUit(sep)}. The situation illus- 
trated in FIG. 7 is a special case of the situation illustrated in 
FIG. 8. 

35 An analysis incorporating the MIT restriction(s) has been 
presented by Grabbe et al in “Modeling and Evaluation of 
Miles-in Trail Restrictions in the National Air Space” 
(A.I.A.A. paper 2003-5628), at the A.I.A.A. Guidance, Navi- 
gation and Control Conference, 11-14 Aug. 2003, Austin, 
Tex., whose content is incorporated by reference herein. In 
40 one embodiment, the analysis models the spacing d ii _ l 
between consecutive aircraft (i and i-1 ) on a route segment as 

4y-=v*-i (dep)-f f _! (dep)), (1 8) 

45 where U(dep) is the actual departure time for aircraft no. k 
(k=i, i-1). This assumes that the time required to reach cruise 
altitude is substantially the same for each of the aircrafts i and 
i-1 and that the true airspeeds for each of the aircrafts i and 
i-1 are substantially the same. Equation (18) can be modified 
to model aircraft separation along a great circle segment, as 

tG_i=(r £ 4-/z i _i) Isin co(r-r 2 -)-sin co(r— r z -_ x ) I , (19) 

® = Vf- 1 /(''£+*i- 1 )» ( 20 ) 

55 where x E is a representative radius of the Earth and h i _ l (=hj 
is the cruise altitude of each aircraft. An analytical miles-in- 
trail (or minutes-in-trail) model works with a MIT time dif- 
ference 

ATL-i =^(dep)-r i _ 1 (dep )=d i /v i _ l , (21) 

60 

and requires that 

ATi^dm/v^, (22) 

where AL is the corresponding MIT minimum separation 
65 distance. This analysis can be extended from two consecutive 
aircraft to N consecutive aircraft (N^2), all traveling the 
same route segment. 
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A second approach for MIT analysis uses a linear program- 
ming model and seeks to minimize a sum 


( N (slots) N(aircraft) (23) 5 

A = mini ^ /^-(slot) - tj(dep)}, 

( i = 1 7=1 


subject to the constraints in Eqs. (22), where N(slots) and 10 
N(aircraft) are the number of aircraft loading slots and the 
number of aircraft, respectively, and nij is a positive weight- 
ing factor (optionally uniform). The weighting factors are 
subject to the following constraints: 


N(slots) 

£ n U = 

1 = 1 

N (aircraft) 

Z "W = 1 

7=1 


In another situation, an aircraft, either enroute or not yet 
departed, inquires about availability of a gate during a 25 
selected time interval, including its estimated arrival time at 
the aircraft’s intended destination. If a landing slot is likely to 
be available for the selected time interval at the selected 
destination, the system advises that the flight can proceed as 
planned. If a landing slot is not likely to be available in the 3Q 
selected time interval at the selected destination, the system 
proceeds as discussed in Section I. 

K. Conflict Detection and Resolution 

FIG. 8 illustrates a spatial relationship between first and 
second aircraft (n=l and n=2) traveling along individual 35 
routes in the same region. Beginning at an initial reference 
location, r=r Q (n=l, 2), and an initial velocity, v=v 0 (n=l, 

2), for each of the aircraft at the same time, t=tO, along the 
respective flight routes, the separation distance 

^(0 = l?*o,i+Vo,i(^-?0)-^- v o,2(^)l (26) 40 


15 


(24) 


(25) 


20 


is computed and minimized with respect to time to determine 
a projected minimum separation distance D(min) given by 

J D(min) 2 ={Ar 1 ; 2 2 Av 1 ? 2 2 -(Ar l 52 -Av l 52 ) 2 }/(Av 1 ? 2 ) 2 , (27) 

45 

Ar lj2 =(Ar 0>1 cos xl cos cos x2 cos X2,r 05 i cos 

xl sin Xl-r 052 cos %2 sin A2,r 0 : sinTl-r 02 sin 
■*2), ’ ’ (28) 

Av lj2 =(r 0JL cos al cos |31-v 02 cos a2 cos |32 ,r 0 l cos 
al sin pi— v 0 2 cos a2 sin p2,v 0 1 sin al-v 0j2 sin 
a2), ’ ’ (29) 

The computed minimum separation time, 

7(min)-ZO=-(Ar i, 2 - Av 12 )/(Av 12 ) 2 , (30) 55 

is required to be non-negative, or the minimum separation 
distance is ignored. 

This minimum separation distance is compared with a 
selected threshold separation distance D(tlir) (typically 3-5 
miles in horizontal separation and 1 000-2000 feet in vertical 60 
separation) to determine if, based upon the projected location 
vectors, the two aircraft will pass too close to each other (i.e., 
D(min)<D(thr)). If the answer to this query is “yes,” one or 
both of these aircraft is advised to alter one or more param- 
eters of its present velocity vector by a selected amount in 65 
order to avoid a separation “incident,” corresponding to 
D(min)^D(thr). If the answer to this query is “no,” the two 


aircraft are allowed to continue, using the present parameter 
values for their velocity vectors. When one or both of the 
aircraft changes at least one velocity vector parameter, either 
sua sponte or in response to a request by the system, a new 
value of D(min) is computed, using the now -modified values 
of the velocity vector parameters, and the comparison process 
is repeated. 

A minimum separation distance D(min) can also be esti- 
mated, using a quadratic or parabolic extension model, rather 
than the linear extension model used in Eq. (26). A flight 
segment of each aircraft is assumed to lie in a plane and to 
approximate a great circle (GC) route, and the location of the 
aircraft is approximated by a quadratic function of the time 
variable t, 


r(Z;app)= \r(t=tO) 1 {ul +a v (Z-zO)+a a (Z-zO) 2 /2} , 

(31) 

a=a vp +a vs , =ul -a vp +u2-a vs , 

(32) 

^a =a ap +aa -as‘ =ul '^p+U^a^, 

(33) 


where ul and u2 are unit length vectors parallel to r(t=t0) and 
to v(t=t0) in the plane GC, respectively, and perpendicular to 
each other. 

The great circle flight route is described by the vector 
equation 

r(t;GC)=\r(t=t0)\{ul cos [<D(f-rf))+<|>]+ 

u2 sin o)([(f-*0)-H|>]} (34) 

where wHv(t=t0)l/lr(t=t0)l and 4> is a phase angle defining an 
initial aircraft location. In the most general case, the vector 
coefficients a vp7 a vs , a ap and a as are determined by minimiz- 
ing an error integral e(t0;T) based on the difference lr(t;app)- 
r(t;GC)l 2 , given by 


£(/0; T) = 


(35) 


f 


|r(r = r0)| 2 |wl-{l -cosw(r- tO) + a vp {t-tO) + a ap (t- tO) 2 / 2} + 


{w2{-sinw(r - /0) + or vs (r - *0) + a as (t - zO ) 2 /2}} 2 dt 


Taking account of the perpendicularity of the vectors ul and 
u2, the minimization equations become 


dsjda vp = I" |r(f = f0)| 2 {l - cos[w(r- rO) + <p] + 

JtO 


(36A) 


9 (%ap — I 

JtO 


ds/da ap = I |r(r = /0)| 2 {1 - cos[w(r- rO) + <p] + 


(36B) 


U'vpK 1 IU / “ 


dsj da vs = (36C) 

f |r(r = r0)| 2 {-sin[w(r -t0) + <p] + 2 a vs (t - tO) + a as (t - tO) 2 / 2} 

JtO 

= 0 


dsj da as = 

f |r(r = f0)| 2 {-sin[w(r - *0) + <f>] + a vs (t - zO) + 2 a as (t - tO) 2 / 2} 
JtO 

(t- tO) 2 /2dt = 0 


(36D) 
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Equations (36A)-(36D) provide two pairs of coupled equa- 
tions: 


A t 

B i 

(X V p 

= 1 Cl 1 

| A2 

B2\ 

a ap 

= | C2 | 

1 A3 

B3\ 

1 ®vs 1 

= CJ 

| A4 

B4\ 

1 Q-as I 

= I c4 r 


A 1 =$to T \ r(t=t§) l 2 {2 (7-*0) 2 }dt, 

A2=J i0 T \r(t=iO) \ 2 {2(t-t0) 2 /2}dt, 

A3=f tO r \r(t=t0) \ 2 {2(t-tO) 2 Jdt, 

A4=f i0 T \r(t=tO) \ 2 {(t-lQf !2}dt, 

Bl=$J\r{t=tQ>)\ 2 {{t-1f)?!2}dt, 

B2=f tO r \r(t=t0) \ 2 {(t-tfff}dt, 

B3=Sj\r{t=tO)\ 2 {{t-lQfl2}dt, 

B4=$J\r{t=t0)\ 2 {{t-tnf}dt, 

Cl=/ /0 r lr(?=rt))l 2 { 1-cos [o)(f-ZO)-H|>]}(f-/0)rfr, 

C2=/ /O r lr(f=/0)l 2 { 1-cos [to{t-ttt)+ty\{t-tti) 2 dt/2, 
C3=J^ r lr(r=/0)l 2 {-sin [<D(f-tf)4<|>](f-fO)rff, 

C4=h 0 r lr(f=tf))l 2 {-sin [to(t-1tt)+<k\{t-t()) 2 dt/2. (37C) 

The minimum separation distance D(min) for two aircraft 
(numbered k=l, 2), whose location vectors are approximated 
as in Eq. (3 1 ), is determined by solving a cubic equation in the 
variable t-tO, namely 

2Ar--Av+2{AvAv+2ArAfl}(if-rt))+6AvAab-/0) 2 + 

4Aa-Av(r-A)) 3 =0, (38) 

where Ar, Av and Aa are the vector differences for the location 
r, velocity v and acceleration a for the two aircraft at t=tO, 
determined using Eqs. (31)-(33). Several straightforward and 
simple methods are available for solving cubic equations, 
such as Eq. (38). A numerical solution (t-tO=t so/ ) is inserted 
into an error term 

e(min)=l Ar+Av-f soi +Aa it sol ) 2 1 2 , (39) 

and this error term is compared with a threshold value D(thr) 2 
to determine if a conflict of the two aircraft is predicted to 
occur. This great circle approximation can also be used for 
trajectory prediction. 

K. D. Bilimoria, in “A Geometric Optimization Approach 
to Aircraft Conflict Resolution” (A.I.A.A. Paper 2000-4265), 
A.I.A.A. Guidance, Navigation and Control Conference, 
14-17 Aug. 2000, Denver, Colo., sets forth an optimized 
method for resolution of an aircraft “conflict,” defined as a 
situation in which two aircraft moving in a common (hori- 
zontal) plane, are projected to pass within a threshold distance 
D(thr) of each other. The content of this article is incorporated 
by reference herein. Conflict detection may use linear or 
nonlinear trajectory prediction. Given two aircraft, A and B, 
spaced apart by a distance r LOS , and a velocity v rel of A 
relative to B, a conflict is predicted to occur if the predicted 
relative trajectory of A (A moving relative to B) will pass 
through at least one point of a sphere S(B), or circle in two 
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dimensions, centered at B and having a radius D(thr), as 
illustrated in FIG. 9. This conflict condition is expressed as 

Dimm^rzoslsmiXLos-Xrei) K^(thr), 

(40) 

fLOS={(XB-XA) 2 +(y3-yA) 2 } 1/2 , 

(41) 

Kei\={vA 2 W2v A v B COS(x 4 -X,)} 1/2 

(42) 

XLas^tair 1 { (y B -y A )l{x B -x A )} 

(43) 

^Urtan- 1 {(v^ sin Xa _ Vb sin Xb )!{(v a cos Xa- 
Xb )}; 

-v B cos 

(44) 


This conflict can be avoided by (1) changing the relative 
heading angle % re of A relative to B to a modified value 

X*rei=XLOs ±s ' m ~ 1 r LOS } , (45) 

corresponding to the relative trajectory of A being tangent to 
the sphere S(B) at one or two surface points, as indicated in 
20 FIG. 9. Where a conflict is present, the relative heading 
change, 

A Xre=X*re~Xre’ ( 46 ) 

is a fundamental parameter, a measure of the change in at least 
25 one trajectory parameter for A and/or B to avoid the predicted 
conflict. 

The conflict can be avoided (1 ) by relative heading change, 
(2) by change of the relative velocity vector \ reh (3) by change 
of a combination of relative heading and relative velocity 
30 vector, (4) by change of altitude of one or both aircraft and/or 
(5) by a change in aircraft ascent rate or descent rate. Where 
relative heading is to be changed, aircraft A and aircraft B can 
be assigned fractional contributions, f A and f B , with f^+f^l , 
to the total relative heading change y* rel , according to a 
35 selected assignment rule. The corresponding fractional 
changes in relative heading become 



XrelA = Xre+fA (X*re~Xre), 

(47 A) 

40 

XrelJ}=Xre+f B (X*re-Xre)- 

(47 B) 


Where a relative heading change is to be made only for 
aircraft A, the corresponding new heading angle is deter- 
mined to be 

45 5C4=X*r^- sin_1 {(V^)sin(x >i, reU -^)} ? (48) 

assuming that the magnitude of the argument of the inverse 
sine function in Eq. (47) is no greater than 1 . 

Where a speed change only is to be implemented, the 
50 modified air speed for aircraft A is determined by 

v^^lsin^rXsysinOc’^rX^)}. (49) 

which is an implicit nonlinear relation between v*^, v B , % A 
and Xs- Equation (49) has two solutions, corresponding to the 
55 two surface tangent points indicated in FIG. 9 . Bilimoria also 
develops an optimal change involving both heading change 
and velocity change. 

L. Direct-to Routing 

60 Direct-to routing is incorporated as an option, to avoid use 
of dog leg route segments between flight route waypoints 1, 2 
and 3, as illustrated in FIG. 10 , when a direct flight from 
waypoint 1 to waypoint 3 is predicted to save at least a 
threshold amount of time At(DTR). Where direct-to routing is 
65 activated, the system estimates the time required for the air- 
craft to travel from waypoint 1 to waypoint 2 to waypoint 3, 
taking account of the local weather, applicable wind field, 
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airspace restrictions and aircraft performance data (“flight 
constraints”). The system then estimates the time required to 
travel from waypoint 1 directly to waypoint 3 (the direct-to 
route), incorporating the corresponding flight constraints and 
compare the estimated times. If the time required to travel the 
conventional route segments (1 to 2 to 3) is at least a selected 
threshold increment At(DTR) (e.g., 60 sec) greater than the 
time required to travel the direct-to route segment (1 to 3), the 
conventional route segments are replaced by the direct-to 
route segment. Otherwise, the flight continues along the con- 
ventional route segments. For each three consecutive way- 
points, this process is optionally repeated. Direct-to routing is 
discussed in H. Erzberger et al, Direct-To Tool for En route 
Controllers,” Proc. IEE Workshop on Advanced Technologies 
and their Impact on Air Traffic Management in the 2\ st Cen- 
tury,” Capri, Italy, 26-30 Sep. 1999 and in B. Sridhar et al, in 
“Benefits of Direct-To Tool in National Airspace System,” 
I.E.E.E. Trans, on Intelligent Transportation Systems, vol. 1 
(2000). The content of these references is incorporated by 
reference herein. The Sridhar et al article applies the 
Erzberger et al model to a particular CTAS site (Fort Worth 
ARTCC), and subsequently to all ARTCC in the NAS, reap- 
plies a modified direst-to routing procedure that is not as 
complex as the CTAS model, and compares the results with 
the corresponding CTAS results. The two models agree 
closely. The modified direct-to routing procedure is part of 
the system disclosed here. 

M. Playbook and CDR Route Evaluation Tools 

The F.A.A. has put together, and continues to revise, a set 
of National Playbook Routes (NPRs), including specified 
waypoints, for a flight between any two of a major East Coast 
airport, a maj or Midwest airport, a maj or Southern airport and 
a major West Coast airport. FIG. 11 illustrates a sequence of 
waypoints between several West Coast airports (LAX, SFO, 
SEA, etc.) and several East Coast airports (JFK, BOS, etc.). 
An NPR route can be specified in a flight plan and used when 
severe weather does not permit a more direct flight by another 
route. For example, a flight from Seattle to Boston that must 
avoid severe weather across the North Central Plains might 
use an NPR route illustrated in FIG. 11. 

Another series of flight routes between a source or origin 
airport and a destination airport is provided by the F.A.A.’s 
Coded Departure Routes (CDRs), provided by the Air Traffic 
Control System Command Center as a sequence of waypoints 
between the source and destination. An example of a CDR 
route between JFK Airport and O’ Hare Airport is shown in 
Table 2 . The CDRs may cover a larger number of airports than 
does the NPR system, and each ARTCC that is traversed by a 
CDR flight route is indicated in this Table. 

The invention allows (1) addition of an aircraft on an NPR 
or CDR and (2) analysis and prediction of NAS-wide impact 
of use of such a route. 

N. System- Wide Optimization 

The system-wide optimization capabilities of the invention 
can be used to calculate an optimal combination of restric- 
tions (i.e. miles -in- trail, minute-in-trail, reroutes, ground 
delay programs and ground stops), which minimize airline 
delays while ensuring that the capacity of scarce NAS 
resources, such as sectors, airports and airways, is met. To 
accomplish this task, detailed models of each of the afore- 
mentioned restrictions are implemented in the invention, for 
example, in connection with miles-in-trail (or minutes-in- 
trail) and rerouting capabilities of the system. The system- 
wide optimization capability can be used in either a “what-if ’ 
mode or a “simulation” mode to perform both real-time plan- 
ning or post-operations analysis studies. 
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In calculating the optimal combinations of restrictions to 
impose, applicable constraints are included to ensure that all 
solutions are equitable from the perspectives of the air carrier 
and the air traffic service provider. In a first example, when 
5 rerouting east-bound traffic around a convective weather cell , 
illustrated in FIG. 12, the invention ensures that traffic is 
equally distributed between the two available routes, labeled 
1 and 2, to ensure that the underlying sectors are not con- 
gested. At the same time, the invention also ensures that no 
10 single airline is forced to fly predominantly along the longer 
and less optimal of the two available routes. 

A second example of the system-wide optimization capa- 
bilities of the invention is illustrated in FIG. 13, where the 
simulation capabilities are used to calculate the NAS -wide 
1 5 impact of varying the departure rates from La Guardia Airport 
(LGA) and Newark Liberty International Airport (EWR) to 
other airports. Because the LGA and EWR airports are adja- 
cent to each other, the cumulative enroute time delays for 
these two airports are not independent of each other. The 
20 dashed line FIG. 13 represents a boundary between those 
airport departure rates that lead to NAS congestion and those 
departure rates that do not. Based on the results presented in 
FIG. 13, the optimal departure rates from LGA and EWR are 
20 and 21 (departures per hour), respectively. This combina- 
25 tion of departure rates ensures that NAS -wide congestion is 
avoided or minimized, while limiting the cumulative airline 
delay to a maximum of 6000 sec. Similar results can be 
generated looking at any combination of restrictions that 
routinely impact congestion and other effects on the NAS. 

30 

O. Overall Procedure 

FIGS. 14 a, 146, 14c and 14 d illustrate a procedure for flow 
of information according to an embodiment of the invention. 
FIGS. 14a and 146 describe the flow of information from air 
35 traffic service provider’s decision-making, and FIGS. 14c and 
14 d describe the flow of information from air traffic service 
user’s decision making. The system first determines, in step 
141, for a given flight or given group of flights, whether the 
flight(s) is active and has a current track and a flight plan or is 
40 based upon a proposed flight plan, which is expected to 
become active at a future time. These data consisting of 
tracks, active flight plans and proposed flight plans are 
recorded, in step 143, and stored in the recorded flight data- 
base (RFDB), in step 145, for use at a later date. Real-time 
45 data from step 141 or historical data from the RFDB are used 
for further processing. The user selects (i) live mode or (ii) 
simulation mode or (iii) playback mode for the flight(s), as 
defined in step 147. In step 149, the system determines if the 
user has selected playback mode. Because only recorded data 
50 can be played back, the playback mode uses data from RFDB . 

If the answer to the query in step 149 is “no,” in step 151 the 
system moves along path 1 and determines, in step 151, if this 
flight(s) is impacted by NAS constraints including one or 
more of the following constraints: playbook routes; GS/GDP 
55 constraints; MIT constraints; local re-routing constraints ; (re) 
sectorization constraints; and departure restrictions. In step 
152, the system allows modification of one or more NAS 
constraints provided in step 151. The system also moves 
along path 5 and provides real-time flight data from step 141 
60 or recorded flight data from RFDB (step 145) to step 182 to 
enable decision-making from air traffic service user’s per- 
spective (discussed in the following). 

One or more defining flight parameters (flight route; depar- 
ture time; flight altitude; flight speed; flight heading; and 
65 destination airport) are modified in step 153 to comply with 
the NAS constraints in step 151. These defining flight param- 
eters are also altered via path 6, as discussed in the following, 
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based on the outcome of collaborative decision-making 
between the air traffic service provider and the air traffic 
service user in step 181 (FIG. 14c). The system then moves 
via path 1 to step 155 to predict flight trajectories (locations at 
future times) of both active aircraft and proposed aircraft, 
using flight parameters from step 153, rapid update cycle 
(RUC) wind velocity forecast data (step 157) and information 
from an aircraft performance database (step 159) containing 
nominal performance data for different types of aircraft. The 
system uses the predicted trajectories to forecast the demand 
for airspace and airport resources, in step 161, where one or 
more of the following quantitative measures of flight activity 
are estimated: traffic count in one or more selected sectors 
(sector count) ; traffic count over one or more fixes (fix count); 
arrival counts at selected airports; departure counts at a 
selected airports; FCA traffic counts; and/or special use air- 
space traffic counts for selected SUAs. Step 161 relies on 
geometric information from an airspace adaptation database, 
provided in step 162. 

If the answer to the query in step 149 is “yes” so that 
playback mode is desired, the system obtains relevant trajec- 
tory information directly from the RFDB (step 145) and fol- 
lows path 2, circumventing the trajectory prediction step in 
155, to forecast demand (step 161). 

Irrespective of the answer to the query in step 149, the 
system then moves to step 163, where a graphical user inter- 
face (GUI) and visualization tools module provide relevant, 
visually perceptible illustrations of aircraft location, flight 
route, severe weather data (step 165), computed demand esti- 
mates (step 161) and demand estimates from an historical 
database (step 167). The system then determines, instep 169, 
if a playback mode was requested earlier in step 149. If the 
answer to the query in step 149 is “yes,” playback is provided, 
based on the presently assembled information, and no further 
action is required (step 171). 

If the answer to the query in step 169 is “no” so that a live 
mode or simulation mode is specified, the system moves to 
step 173 and determines if additional NAS constraints are 
needed for mitigating imbalances between demand for, and 
the available capacities of, the airspace and airport resources, 
in order to manage air traffic. If the answer to the query in step 
173 is “no,” the system applies a conflict detection and reso- 
lution (CD&R) analysis and response to the active and pro- 
posed flights, in step 175, and determines, in step 177, 
whether the flights are conflict-free after application of the 
CD&R analysis and response. 

If the answer to the query in step 173 is “yes,” the system 
follows path 4 and determines one or more of the NAS con- 
straints that need modification (step 152), changes the NAS 
constraints accordingly in step 151, determines which flights 
are impacted by these new NAS constraints in step 151, 
changes one or more of the selected route parameters to 
comply with the new constraints (step 153), and continues 
along path 1 as before. 

If the answer to the query in step 177 is “no,” the system 
moves along path 3 to step 153 and modifies at least one of the 
following flight parameters: flight route; departure time; 
flight speed; altitude; flight heading; and destination airport. 
After step 153, the system again proceeds along path 1. 

If the answer to the query in step 177 is “yes,” the system 
follows path 7 and generates NAS decision data from the 
service provider’s perspective (optionally including a new set 
of NAS constraints and flight parameter changes), in step 179. 
The system continues along path 7 to step 181, where col- 
laborative decision-making between the air traffic service 
provider and the air traffic service user occurs. The system 
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proceeds along path 6 to steps 152 and 153, depending upon 
the results of collaborative decision-making and proceeds 
again along path 1 . 

Service providers such as the Federal Aviation Adminis- 
5 tration (FAA) in the United States would typically perform 
the procedures in steps 141 through 179 in FIGS. 14a-14b. 
The users of air traffic services are typically commercial 
aviation, business aviation, general aviation, military and 
individual pilots. Both air traffic service providers and air 
l o traffic service users (collectively referred to as “users” herein) 
can use the system. 

Along path 7, the system proceeds to step 181, collabora- 
tive decision making and, in parallel, to step 182, where it is 
determined if the air traffic service user’ s flights are impacted 
15 by NAS constraints. Step 182 uses real-time data from step 
141 or historical data from step 145, received via path 5. 
Desired modifications to NAS constraints in step 211 (FIG. 
14 d) are also received in step 182 via path 10. Step 182 is 
substantially similar to step 151. 

20 One or more trajectory alternatives are generated in step 
183, including wind optimal routes and NPR routes and user- 
preferred routes to mitigate the impact of NAS constraints on 
user’s flights. The alternative trajectory generation step 183 
utilizes RUC wind data (step 185) and aircraft performance 
25 data (step 187) that is generic (as in step 159) or is specific to 
user’s particular fleet of aircraft. 

Flight parameters including flight route; departure time; 
flight altitude; flight speed; flight heading; and destination 
airport are modified in step 184 to comply with the proposed 
30 NAS constraints provided in step 182 and to realize the alter- 
native trajectories generated via step 183. Trajectories of both 
active and proposed aircraft are predicted in step 188 using 
the flight parameters specified in step 184, RUC wind velocity 
forecast (step 185) and aircraft performance data (step 187). 
35 The collaborative decision making step often involves 
negotiation between the service provider and the service user 
concerning modification of NAS constraints (step 152) and 
the resulting defining flight parameters (step 153). If, as a 
result of such negotiation, one or more NAS constraints and/ 
40 or one or more defining flight parameters are changed, the 
procedures of steps 151 through 179 are repeated. 

From step 188, the system moves to step 189, demand 
forecasting using aircraft adaptation data (step 190), where 
one or more of the following quantitative measures of flight 
45 activity are estimated: traffic count in one or more selected 
sectors (sector count); traffic count over one or more fixes (fix 
count); arrival counts at selected airports; departure counts at 
a selected airports; FCA traffic counts; and/or special use 
airspace traffic counts for selected SUAs. The procedures in 
50 steps 161 and 189 are substantially identical 

The system then moves to step 191, where a graphical user 
interface and visualization tools module provides relevant, 
visually perceptible illustrations of aircraft location, flight 
route, severe weather data from step 193, computed demand 
55 estimates from step 189 and/or historical airspace demand 
data from database in step 195. The procedures in steps 163 
and step 191 may be substantially the same, or step 191 may 
include additional illustrations especially tailored from the 
airspace service user’s perspective. 

60 The system then moves along path 8 in the following man- 

ner: (1) to step 201 and determines if one or more flights need 
additional modification; and (in parallel) (2) to step 203 and 
determines if one or more of the NAS constraints need addi- 
tional modification. If the answer to the query in step 201 is 
65 “no” so that no additional modifications are needed), the 
system generates user decision data, in step 209, which may 
include proposals for changes in defining flight parameters 
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(step 181). If the answer to the query in step 201 is “yes,” the 
system implements one or more of the following actions, in 
step 207: modify flight route; modify flight departure time; 
cancel a flight; and provide a substitute flight in lieu of the 
cancelled flight. These changes are provided to step 184 via 5 
path 11 for reassessment via modules 184, 188, 189 and 191. 

If the answer to the query in step 203 is “no,” the system 
moves to step 209 to generate and present user decision data, 
which may include proposals for changes in NAS constraints io 
(step 181). If the answer to the query in step 203 is “yes,” the 
system proposes modifications in one or more NAS con- 
straints, in step 211, and provides these data to module 182 
via path 10. The impact of the proposed modifications to the 
NAS constraints can be reexamined via modules 182, 183, 15 
184, 188, 189 and 191 along with the supporting data modules 
185, 187, 190, 193 and 195. Once the desired set of proposed 
NAS constraints and flight parameters is obtained by repeated 
reevaluation via paths 11 and 10, the system then moves to 
step 209, then to step 181, where both the service provider and 20 
the service user, or several users, collectively agree on the 
choice of NAS constraints and flight parameters. These 
agreed upon choices are then realized in steps 152 and 153. 
The procedures illustrated in FIGS. \4a-\4d are applied to 25 
one or more aircraft flights and to the corresponding aircraft. 

The overall system-procedure, illustrated in one embodi- 
ment in FIG. 14, may use information and features from the 
graphical user interface (GUI), the weather and winds data 
module, the weather/winds interpolation module, the filed 30 
flight plans module, the aircraft performance database, the air 
traffic monitoring module, the route parser and/or trajectory 
predictor module, the traffic analyzer module, the miles-in- 
trail and/or minutes-in trail restriction module, the conflict 
detection and resolution (CD&R) module, the direct-to mod- 35 
ule, the playback and CD&R evaluation module, and/or the 
system-wide optimization module, as discussed in the pre- 
ceding Sections, A, B, C, D, E, F, G, H, I, J, K, L, M and N. 

40 

TABLE 1 

Aircraft Performance Data. 


CC B757 PERFORMANCE FILE Oct. 1, 1998 

CC 

CC AC/Type: B757. DADS Revision: 1 . 1 

CC 

CC Units: 

CC Speeds: CAS Mach Mass (kg) Temperature: ISA 

CC 

cas_climb = 290; 

mach climb = 0.78; 

CC 

cas_cruise = 320; 

mach cruise = 0.80; 

CC 

cas_descent = 300; 
mach_descent = 0.78; 

CC 

mass_low = 69600; 
mass_fain = 95000; 
mass_high = 110000; 

CC 

max_alt =42000; 

CC 

cruise_alt = 37000; 
cruise_alt_east = 37000; 
cruise_alt_west = 35000; 

CC 

CC cruise data = [FL; TAS (knots); fuel low(kg/min); nom(kg/min); 
high(kg/min)] 


45 


50 


55 



TABLE 1 -continued 

Aircraft Performance Data, 
fuel consump. 

FL TAS (low) (med) (high) 

30 261 40.8 51.2 58.9 

40 265 40.9 51.4 59.1 

60 272 41.1 51.7 59.5 

80 280 41.3 52.0 59.8 

100 289 41.5 52.3 60.2 

120 379 58.3 65.5 70.7 

140 390 58.5 65.8 71.1 

160 401 58.7 66.1 71.5 

180 413 58.9 66.4 71.9 

200 425 59.0 E6.7 72.3 

220 438 59.2 66.9 72.7 

24 451 59.3 67.2 73.1 

260 465 59.4 67.5 73.5 

280 475 58.7 67.0 73.2 

300 471 55.2 64.3 71.1 

320 467 52.2 62.1 69.5 

340 463 49.5 60.5 68.5 

360 459 47.3 59.3 68.1 

380 458 45.7 58.9 68.6 

400 458 44.5 58.9 69.6 

420 458 43.3 58.9 70.6 

CC 

CC climb data = [FL; TAS(knots); ROCD low(fpm); nom(fpm); 

high(fpm); fuel nom(kg/min)] 


FLSR 

TAS 

ROCD(low) 

(med) 

(high) 

fuel consump. 
(med) 

0 

169 

3760 

2320 

1730 

170.4 

5 

170 

3740 

2300 

1710 

169.1 

10 

171 

3720 

2280 

1690 

167.7 

15 

172 

3700 

2270 

1670 

166.3 

20 

174 

3680 

2250 

1660 

164.9 

30 

261 

5120 

3460 

2800 

173.9 

40 

265 

5060 

3410 

2750 

171.1 

60 

272 

4910 

3290 

2650 

165.5 

80 

280 

4770 

3180 

2540 

159.9 

100 

289 

4610 

3050 

2430 

154.4 

120 

344 

4670 

3150 

2540 

154.2 

140 

354 

4470 

2990 

2400 

148.6 

160 

365 

4260 

2820 

2250 

143.0 

180 

376 

4030 

2650 

2100 

137.4 

200 

387 

3800 

2480 

1940 

131.7 

220 

399 

3570 

2300 

1780 

126.1 

240 

412 

3320 

2110 

1610 

120.5 

260 

425 

3070 

1910 

1440 

114.8 

280 

438 

2810 

1710 

1260 

109.1 

300 

452 

2540 

1510 

1080 

103.4 

320 

455 

3210 

1830 

1240 

97.0 

340 

451 

2880 

1540 

970 

90.1 

368 

447 

2540 

1230 

670 

83.4 

380 

447 

2010 

850 

330 

76.9 

400 

447 

1680 

540 

30 

70.5 

420 

447 

1350 

230 

0 

64.1 

CC descent 
fuel 

_data = 

[FI; TAS (knots); ROCD(fpm); ftiel (kg/min)] 

FLSK 

TAS 

ROCD 

consumed 



0 

132 

1340 

19.0 



5 

133 

1350 

18.8 



10 

134 

1360 

18.7 



15 

151 

1280 

18.5 



20 

193 

1210 

18.3 



30 

217 

1250 

18.0 



40 

241 

1340 

17.7 



60 

272 

1490 

17.0 



80 

280 

1520 

16.2 



100 

289 

1560 

15.7 



120 

356 

2020 

15.0 



140 

366 

2060 

14.3 



160 

377 

2090 

13.7 



180 

388 

2120 

13.0 



20.0 

400 

2160 

12.3 




65 
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TABLE 1 -continued 


Aircraft Performance Data. 



3 JFKORD64 KJFK 


4 JFKORD80 KJFK 


5 JFKORD95 KJFK 


6 JFKORDCA KJFK 


7 JFKORDDJ KJFK 


8 JFKORDJ6 KJFK 
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TABLE 2 -continued 


Route Departure Route Departure Arrival Traverted 

# Code Origin Destination Fix String ARTCC ARTCC ARTCCs 


9 JFKORDJV KJFK KORD 


10 JFKORDP5 KJFK KORD 


1 1 JFKORDPH KJFK KORD 


1 2 JFKORDRF KJFK KORD 


13 JKFORDX6 KJFK KORD 


14 JFKORDXU KJFK KORD 


J6 

COLNS 

J134 

FLM 

J24 

VHP 

OKK 

OKK1 

KORD 


GREKI 

KJFK 

ZNY 

ZAU 

CZY 


GREKI 



ZAU 


V419 



ZBW 


JUDDS 



ZMP 


CAM 

ART 

YCF 

YEE 

ASP 

TVC 

GRB 

MSN 

JVL 

JVL4 

KORD 



ZNY 

RBV 

KJFK 

ZNY 

ZAU 

ZAU 


RBV 



ZID 


J230 



ZKC 


AIR J80 



ZNY 


CAP 

PNT 

V227 

PLANO 

KORD 



ZOB 

COATE 

KJFK 

ZNY 

ZAU 

ZAU 


COATE 



ZNY 


J36 

FNT 

PMM4 

KORD 



ZOB 

WAVEY 

KJFK 

ZNY 

ZAU 

ZAU 


WAVEY 



ZBW 


EMJAY 



ZDC 


J174 



ZID 


ORF 



ZNY 


PSKIIU 

VHP 

OKK 

OKK1 

KORD 



ZTL 

RBV 

KJFK 

ZNY 

ZAU 

ZAU 


RBV 



ZDC 


J230 



ZID 


SAMME 

J6 

EYTEE 

J149 

FWA 

DXI3 

KORD 



ZNY 

GREKI 

KJFK 

ZNY 

ZAU 

CZY 


GREKI 



ZAU 


V419 



ZBW 


JUDDS 



ZNY 


CAM 



ZOB 


J547 

SYR 

J63 

EHMAN 

YXU 

J547 

PMM 

PMM4 

KORD 
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What is claimed is: 

1. A method for estimating a minimum distance of 
approach of two aircraft that are airborne, the method com- 
prising: 

providing information on an initial location vector rO(t=tl ; 
n), on an initial velocity vector rl(t=tl;n) and on an 
initial acceleration vector r2(t=tl ;n) at a selected time 
t=t 1 , for each of N aircraft, numbered n=l , . . . , N (Ni^2) 
that are airborne; 

approximating said location vector r(t;n) for aircrafts num- 
ber n=nl and n=n2 (nl ^n2) over a selected time interval 
[tl,t2] by quadratic vector functions of time, 

r(t;n 1 ;app)=rO(« l)+rl (i n 1 )■(/-/! )+r2{nl)-{t-tl ) 2 , 

r(t;n2;app)=r0(n2)+rl(n2y(t-tl)+r2(n2y(t-tl) 2 , 

Ar(t;app)=r(t;n 1 ;app)-r(U«2;app)=AK)+Arl (r-rl)+ 

Ar2(t-t\f, 

respectively, where 1 1 is a selected time within a selected time 
interval [T1,T2], each of the location vectors r(t;nl;app) and 
r(t;n2;app) substantially describes motion on a great circle in 
a plane, and the vector coefficients rO;nl), rl(nl), r2(nl), 
r0(n2), rl(n2) and r2(n2) are chosen to optimally match the 
vector functions r(t;nl;app) and r(t;n2;app) in the selected 
time interval [T1,T2]; and 

estimating a minimum distance of approach d(min) for a 
magnitude Ir(t;nl)-r(t;n2)l of a vector difference, by 
identifying at least one real time t(min) for which a time 
derivative of the quantity Ir(t;nl)-r(t;n2)l 2 is zero, 

2ArOArl+{Arl-Arl+2ArO-Ar2)(t-tl)+6Arl-Ar2(t-tl) 2 + 

4Ar2 ■ Ar2 (t-t\ ) 3 =0, 

and by interpreting the vector magnitude lr(t=t(min);nl)-r 
(t=t(min);n2)l as the minimum distance d(min). 

2. A method for managing aircraft traffic, the method com- 
prising: 

providing information on location vector r w (t=tm) and 
velocity vector v w (t=tm) for each of N aircraft, num- 
bered n=l, . . . , N (Ni=2) that are airborne and are 
located within a selected air route traffic control center 
(ARTCC), for at least one selected time t=tm, where 
each of the N aircraft is assigned to at least one ARTCC 
sector, numbered s=l, . . . , S (S^2) in the selected 
ARTCC; 

at a time, t=tm'>tm, altering at least one boundary of each 
of at least two selected adjacent ARTCC sectors, num- 
bered s=sl and s=s2 (sl*s2), within the selected 
ARTCC to provide altered sectors, numbered s=sl' and 
s— s2', respectively, where the union of the at least two 
selected adjacent sectors encloses the union of the at 
least two altered sectors; and 

providing information on location vector r„(t=tm’) and 
velocity vector v w (t=tm') for each of the N aircraft, that is 
airborne and is located within the selected ARTCC, for 
the time t=tm’, where each of the N aircraft is assigned to 
at least one ARTCC sector, numbered s=l, . . . , S (S^2) 
in the selected ARTCC. 

3. A method for managing aircraft traffic, the method com- 
prising: 

providing information on an initial location vector r 0w 
(t=tO) and an initial velocity vector v 0w (t=tO) for each of 
N aircraft, numbered n=l, . . . , N (N^2) that are air- 
borne; 

approximating the location vector r(t;n) for the aircraft 
number n=nl and for the aircraft n=n2 (nl^n2; nl, 
n2^N) by vector functions that are at least quadratic in 
a time variable t, 

r(t;n 1 ;app )=rO(t;n 1 )+rl ( t;n 1 ) ■ ( t-t0)+r2 ( t;n 1 )(t-tO) 2 , 
r(t;n2;app)=r0(t;n2)+rl(t;n2y(t-t0)+r2(t;n2)(t-t0) 2 , 
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respectively, relative to a selected initial time tO, where each 
of the location vectors r(t;nl) and r(t;n 2 ) substantially 
describes motion on a great circle in a plane, and the vector 
coefficients rO(t;nl), rl(t;nl), r 2 (t;nl), r 0 (t;n 2 ), rl(t;n 2 ) and 
5 r 2 (t;n 2 ) are chosen to optimally match the vector functions 
r(t;nl) and r(t;n 2 ) in a selected time interval [tl ,t 2 ] ; and 

determining a descent-start location at which the aircraft 
should begin its descent toward a destination airport 
along a substantially linear descent path, where an alti- 
1 tude descent rate experienced by the aircraft i s a selected 

fraction f ( 0 <f^i 1 ) of a maximum altitude descent rate. 

4 . A method for managing aircraft traffic, the method com- 
prising: 

providing information on an initial location vector r 0n 
(t=tO) and an initial velocity vector v 0 „(t=tO) for each of 
15 N aircraft, numbered n=l, . . . , N (Ni^2) that are air- 
borne; 

approximating the location vector r(t;n) for the aircraft 
number n=nl and for the aircraft n=n 2 (nl^n 2 ; nl, 
n 2 ^N by vector functions that are quadratic in a time 
20 variable t, 

r(t;n 1 ;app)=rO(/;H l)+rl ( t;n 1 ) ■ (t-t0)+r2 ( t;n ly(t-tff) 2 , 

r(t;n2;app)=r0(t;n2)+rl(t;n2)-(t-t0)+r2(t;n2)(t-t0) 2 , 

25 respectively, relative to a selected initial time tO, where each 
of the location vectors r(t;nl) and r(t;n 2 ) substantially 
describes motion on a great circle in a plane, and the vector 
coefficients rO(t;nl), rl(t;nl), r 2 (t;nl), r 0 (t;n 2 ), rl(t;n 2 ) and 
r 2 (t;n 2 ) are chosen to optimally match the vector functions 
30 r(t;nl) and r(t;n 2 ) in a selected time interval [tl ,t 2 ] ; 

providing an estimate of a wind velocity vector v w =(v w cos 
0 w ,v w sin 0 W ) at a specified location, where v w is an 
estimated magnitude of the wind velocity vector and 0 W 
is an angle of the wind velocity vector measured relative 
to a selected reference line or reference plane; 

35 providing an estimate of a desired angle Q des of travel of 
said aircraft in an environment including the estimated 
wind velocity vector; and 

orienting a velocity vector associated with said aircraft at 
an angle Q comp relative to the reference line or plane, 
40 where 6 comp is determined by 

tan Q comp ={v des sin 6 de -v w sin Q w }/{v des cos Q des - 
v w scos 0 W }. 

5. A method for managing aircraft traffic, the method com- 
45 prising: 

providing information on an initial location vector r 0w 
(t=tO) and an initial velocity vector v 0 w (t=tO) for each of 
N aircraft, numbered n=l, . . . , N (Ni=2) that are air- 
borne; 

50 approximating the location vector r(t;n) for the aircraft 
number n=nl and for the aircraft n=n 2 (nl*n 2 ; nl, 
n2^=N) by vector functions that are quadratic in a time 
variable t, 

r{t;n 1 ;app)=rO(/;K 1 )+r 1 (t;n 1 ) ■ {t-tff)+r2 ( t;n 1 ) ■ (t-tO) 2 , 

55 

r(t;n2;app)=r0(t;n2)+rl(t;n2y(t-t0)+r2(t;n2)(t-t0) 2 , 

respectively, relative to a selected initial time tO, where each 
of the location vectors r(t;nl) and r(t;n 2 ) substantially 
describes motion on a great circle in a plane, and the vector 
60 coefficients rO(t;nl), rl(t;nl), r 2 (t;nl), r 0 (t;n 2 ), rl(t;n 2 ) and 
r 2 (t;n 2 ) are chosen to optimally match the vector functions 
r(t;nl) and r(t;n 2 ) in a selected time interval [tl ,t 2 ] ; 

providing an estimate of a wind velocity vector v w =(v w cos 
0 w ,v w sin 0 W ) at a specified location, where v w is an 
65 estimated magnitude of the wind velocity vector and 0 w 
is an angle of the wind velocity vector measured relative 
to a selected reference line or reference plane; 
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providing an estimate of a desired angle § des of travel and a 
desired magnitude of velocity of travel v des of said air- 
craft; 

providing an estimate of a desired magnitude of velocity of 
travel v des in an environment including the estimated 
wind velocity vector; and 

providing said aircraft with an associated magnitude of 
aircraft velocity v comp , in the absence of said wind veloc- 
ity vector, that is given by 

V comp = { V des 2+V w 2 -2Vjes vw COS(0^-0 W )} 1/2 . 

6. A method for managing aircraft traffic, the method com- 
prising: 

receiving and storing in a database an estimated location 
vector for a sequence of times over a time interval that 
includes at least at two flight days; and 
using information in the database to estimate a number of 
flights in a selected region, including at least one iden- 
tified ARTCC sector, for at least one prediction time that 
is not included in the at least two flight days; 
estimating a number of aircraft that will be located in the 
selected region at each of a second selected sequence of 
times; and 

when the at least one identified ARTCC sector will contain 
more than a selected threshold number of the aircraft at 
an identified time among the second sequence of times, 
changing at least one boundary between the at least one 
identified ARTCC sector and an adjacent ARTCC sector 
to reduce the number of aircraft contained in the at least 
one identified ARTCC sector at a time preceding the 
identified time; and 

displaying a selected area including the selected region, 
after the at least one boundary is changed, in a visually 
distinguishable format, when the selected region will 
contain no more than the selected threshold number of 
the aircraft at the identified time. 

7. A method for managing aircraft traffic, the method com- 
prising: 

receiving and storing in a database an estimated location 
vector and estimated velocity vector for each of N air- 
craft (Ni^2) at each of a selected sequence of times over 
a time interval that includes at least at two flight days; 
using information in the database to estimate a number of 
flights within an identified ARTCC sector, including a 
selected airport at which the N aircraft are expected to 
land, for at least one prediction time that is not included 
in the at least two flight days; 
estimating a number of the N aircraft that will descend and 
land at the selected airport within each of a selected 
sequence of time intervals; and 
estimating a demand on at least one of (i) at least one 
runway at the selected airport and (ii) a selected group of 
arrival -departure gates at the selected airport. 

8 . A method for managing aircraft traffic, the method com- 
prising: 

receiving and storing in a database an estimated location 
vector and estimated velocity vector for each of N air- 
craft (N^2) at each of a selected sequence of times over 
a time interval that includes at least at two flight days; 
using information in the database to estimate a number of 
flights within an identified ARTCC sector, including a 
selected airport at which the N aircraft are initially 
located, for at least one prediction time that is not 
included in the at least two flight days; 
estimating a number of the N aircraft that will take off and 
ascend from the selected airport within each of a 
selected sequence of time intervals; 
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estimating a demand on at least one of (i) at least one 
runway at the selected airport and (ii) a selected group of 
arrival -departure gates at the selected airport; 
providing information on an initial location vector r 0n 
5 (t=tO) and an initial velocity vector v 0 „(t=tO) for each of 
N aircraft, numbered n=l, . . . , N (N^2) that are air- 
borne; 

aircraft at a time t(est) that is displaced from an initial time 
tO by approximate time increments Atm, numbered m= 

; 1, . . . , M (M^2), where 0<Atl<At2^ . . . <AtM; and 

where the magnitude lr 01 — r 02 l of the difference of the loca- 
tion vectors of aircrafts number n=l and n=2 is esti- 
mated to be less than a selected difference value for at 
least one of the time increments Atm, assigning a con- 
flict avoidance response to at least one of the aircrafts 
15 number n=l and n=2 so that, with the conflict avoidance 
response implemented, the magnitude of the difference 
of the location vectors of aircrafts number n=l and n=2 
for each of the time increments Atm is no less than the 
selected difference value. 

20 9. A method for estimating a minimum distance of 

approach of two aircraft that are airborne, the method com- 
prising: 

providing information on an initial location vector rO(n)=r 
(t=tl ;n) and an initial velocity vector vO(n)=v(t=tl ;n) for 
25 each of N aircraft, numbered n=l , . . . , N (Ni=2) that are 

airborne; 

estimating a location separation vector Ar(t;n) for each of at 
least two aircraft, number n=nl and n=n2 (nl^n2), as 
Ar(t;n)=rO(n)+vO(n)(t-t 1 ) for a selected reference time 
tl; 

estimating a time at which a minimum distance of separa- 
tion occurs for the aircraft, n=nl and n=n2, to be about 
At(min)=tl-(Ar 1 2 *Av 9 2)/)/(Av 12 ) 2 , where Ar 12 =rO(nl)- 
r0(n2) and Av 12 =v0(nl)-v0(n2); and 
estimating a minimum distance of separation to be about 

35 IAr(min)l={Ar lj 2 2 Av lj 2 2 -Ar l5 2'Av 1 2)}/(Av l5 2) 2 . 

10. The method of claim 1, further comprising: 

when said minimum distance of approach is less than a 
selected threshold number for a value of said time t=t 
(min) within said selected time interval [tl ,t2], advising 
40 at least one of said aircraft number nl and number n2 
interpreting this condition as indicating that an aircraft 
conflict is likely to occur between said aircrafts n=nl 
and n=n2 within said selected time interval; and 
when an aircraft conflict is determined to be likely to occur, 
45 allowing at least one of said aircraft number n=nl and 
n=n2 to adopt a conflict avoidance response from a 
group of avoidance responses comprising: (1) changing 
a heading angle for said at least one of said aircrafts 
number n=l and n=2; (2) changing a velocity vector for 
said at least one of said aircrafts number n=l and n=2; 
(3) changing a magnitude of at least one of said initial 
velocity vectors v 01 and v 02 ; (4) changing an altitude of 
flight for said at least one of said aircrafts number n=l 
and n=2; and (5) changing at least one of an aircraft 
ascent rate and an aircraft descent rate for at least one of 
55 said aircrafts number n=l and n=2. 

11 . The method of claim 1 , fiirther comprising providing a 
desired route for at least one aircraft that is substantially at 
least one of a wind-optimal route and an NPR route. 

12 . The method of claim 1 , further comprising choosing 
60 said selected difference value to lie in a range of 3-5 nautical 

miles in a horizontal direction and 1 000-2000 feet in a vertical 
direction. 



